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'  The  asymptotic  stress  and  deformation  fields  of  a  crack  propagating  steadily  and 

quasi-stadcally  into  an  elastic-plastic  material,  characterised  by  Jj-  flow  theory  with  linear 

r —  « 

strain-hardening,  were  first  determined  by  Amazigo  and  Hutchinson  (1977)  for  the  cases  of 
mode  IQ  and  mode  I  (plane  strain  and  plane  stress).  Their  solutions  were  approximate  in 

that  they  neglected  the  possibility  of  plastic  reloading  on  the  crack  faces.  This  effect  was 

Y A«-  * .» 

taken  into  account  by  PonteCastaneda  (1985c)>  who  formulated  the  (eigenvalue)  problem 
in  terms  of  a  system  of  first  order  O.D.E.’s  in  the  angular  variations  of  the  stress  and 
velocity  components.  The  strength  of  the  singularity,  serving  as  the  eigenvalue,  and  the 
angular  variations  of  the  fields  were  determined  as  functions  of  the  hardening  parameter. 

The  above  analysis  does  not  determine  the  amplitude  factor  of  the  near-tip 
asymptotic  fields,  or  plastic  stress  intensity  factor.  In  this  work,  a  simple,  approximate 
technique  based  on  a  variational  statement  of  compatibility  is  developed  under  the 
assumption  of  small-scale  yielding.  A  trial  function  for  die  stress  function  of  the  problem, 
that  makes  use  of  the  asymptotic  information  in  the  near-tip  and  far-field  limits,  is 

postulated.  Such  a  trial  function  depends  on  arbitrary  parameters  that  measure  the  intensity 

<\ 

of  the  near-tip  fields  and  other  global  properties  of  the  solution.  Application  of  the 
variational  statement  then  yields  optimal  values  for  these  parameters,  and  in  particular 
determines  the  plastic  stress  intensity  factor,  thus  completing  the  knowledge  of  the  near-tip 
asymptotic  fields.  The  results  obtained  by  this  novel  method  are  compared  to  available 
finite  element  results. 


MW 


1.  Introduction 


In  the  propagation  of  cracks  in  ductile  materials,  a  stage  of  stable  crack  growth, 
where  the  efects  of  inertia  are  unimportant,  is  often  observed  under  continuously  increasing 
external  load,  and  immediately  following  the  onset  of  crack  growth.  This  phenomenon  is 
intrinsic  to  the  plastic  deformation  at  the  tip  of  a  moving  crack:  as  the  crack  propagates  in 
an  elastic-plastic  material,  it  leaves  behind  a  wake  of  residual  plastic  deformation,  which 
typically  accounts  for  all  of  the  energy  available  from  the  external  source,  and  thus  reduces 
the  tendency  of  the  crack  to  grow. 

Under  small-scale  yielding,  which  assumes  that  the  size  of  the  active  plastic  zone 
surrounding  the  crack  tip  is  small  compared  to  the  crack  size  and  the  overall  dimensions  of 
the  specimen  in  consideration,  an  appropriate  measure  of  the  external  loading  is  the  stress 
intensity  factor  of  the  enclosing  elastic  singularity  fields,  hi  accord  with  the  discussion  in 
the  previous  parragraph,  it  has  been  postulated  that  the  graph  of  such  elastic  stress  intensity 
factor  versus  the  amount  of  crack  growth,  called  tire  resistance  curve,  is  a  material  property 
for  small-scale  yielding. 

Central  to  these  notions,  and  to  die  continued  development  of  fracture  mechanics  in 
general,  is  the  knowledge  of  the  stress  and  deformation  fields  near  die  tip  of  the  advancing 
crack,  and  important  advances  have  been  made  in  recent  years. 

Most  of  die  work,  however,  has  been  in  the  context  of  the  ideal  elastic-plastic  solid. 
McClintock  (19S8)  and  McClintock  and  Irwin  (1965)  were  the  first  to  demonstrate  that  the 
key  to  stable  crack  growth  lies  in  the  nature  of  the  plastic  deformation,  and  found  that  the 
strains  in  the  active  plastic  zone  near  the  tip  of  the  crack  depend  logarithmically  on  the  radial 
distance  from  the  crack  tip,  r.  Later,  Chitaley  and  McClintock  (1971)  gave  the  correct 
asymptotic  solution  to  the  problem  of  a  steadily  and  quasi-statically  growing  crack  in  an 
ideal  elastic-plastic  material  with  Mises  yield  criterion.  Rice  (1974)  demonstrated  that  the 
same  logarithmic  dependence  of  the  strains  on  r  would  be  expected  for  the  cases  of  plane 


stress  and  plane  strain.  Slepyan  (1974)  obtained  the  first  correct  asymptotic  solution  of  the 
plane  strain  problem  in  modes  I  aiw?  n  using  a  Tresca  yield  criterion.  Rice  et  al.  (1980), 
following  Rice  and  Sorensen  (1978),  and  Gao  and  Hwang  (1980),  independently,  solved 
the  plane  strain  problem  in  mode  I  using  a  Mises  yield  criterion,  but  the  validity  of  their 
solution  was  restricted  to  the  case  of  v  -  1/2.  Somewhat  later,  Drugan  et  al.  (1982) 
generalised  the  results  to  the  case  of  v  #  1/2.  The  plane-stress  mode  I  problem  has 
remained  most  elusive,  but  recently,  Ponte  Castafieda  (1985a)  produced  a  solution  to  the 
corresponding  mode  II  problem.  In  addition,  a  general  discussion  of  the  subject  can  be 
found  in  Rice  (1982). 

The  corresponding  work  for  the  hardening  elastic-plastic  solid  is  at  a  less  developed 
stage,  and  most  of  this  work  relates  to  the  linear-hardening  case.  Slepyan  (1973),  using  a 
modified  J2-defarmation  theory,  produced  the  first  asymptotic  solution  to  the  problem  of  a 

mode  in  crack  propagating  steadily  in  a  linear-hardening  material.  Much  later,  Ponte 
Castaneda  (1985b)  extended  Slepyan's  work  to  the  mode  I  case.  Amazigo  and  Hutchinson 
(1977),  using  J2-flow  theory,  produced  steady-state  asymptotic  solutions  for  the  moving 

crack  in  anti-plane  strain,  as  well  as  in  plane  strain  and  plane  stress  mode  I.  For  a  linear- 
hardening  material,  with  Young's  modulus,  E,  and  tangent  modulus,  they  found  a 

variable-separable  solution  with  stress  and  strain  fields  having  a  power  singularity  of 
strength,  s,  which  is  a  function  of  the  hardening  ratio  a  =  E/E,  and  vanishes  in  the 

perfectly-plastic  limit.  These  solutions,  however,  neglected  the  possibility  of  plastic 
reloading  on  the  crack  flanks  which  is  known  to  occur  in  the  elastic-perfectly  plastic 
solutions  mentioned  above.  In  a  recent  publication,  Zhang  et  al.  (1983)  included  reloading 
(as  well  as  a  Bauschinger  effect)  in  their  solution  of  the  plane-strain  mode  I  problem.  More 
recently,  but  independently,  Ponte  Castaneda  (1985c)  formulated  the  problem  efficiently  in 
terms  of  a  system  of  first  order  OD.E.'s  in  the  angular  variations  of  the  components  of  the 
stress  and  velocity  fields,  and  included  as  well  the  effect  of  reloading  in  the  solution  of  the 
anti-plane  strain,  plane  strain  and  plane  stress  (modes  I  and  II)  problems. 


These  asymptotic  solutions  leave  undetermined  certain  parameters  in  the 
expressions  for  the  stress  and  deformation  fields  near  the  tip  of  the  moving  crack.  In  the 
perfectly-plastic  case,  the  parameter  in  question  is  die  size  of  the  plastic  zone  ahead  of  the 
crack,  but  this  parameter  does  not  appear  in  the  lowest  order  term  of  the  asymptotic 
expansion  of  the  fields  as  r  — » 0,  and  its  determination  is  not  critical.  On  the  other  hand,  in 
the  linear-hardening  case,  by  analogy  with  the  elastic  problem,  the  amplitude  of  the  lowest 
order  term  of  the  asymptotic  expansion  of  the  stress  and  deformation  fields,  as  measured 
by  the  plastic  stress  intensity  factor,  is  undetermined  from  the  asymptotic  analysis,  and  its 
determination  introduces  an  important  problem. 

In  this  paper  we  use  the  near-tip  information  known  from  the  work  of  Ponte 
Castaneda  (1985c)  to  complete  the  knowledge  of  the  asymptotic  fields  of  a  steadily 
propagating  crack  in  a  linear-hardening  material  under  small-scale  yielding  conditions  by 
determining  the  plastic  stress  intensity  factor  of  these  near-dp  fields  as  a  function  of  the 
hardening  parameter.  Other  relevant  information,  such  as  the  size  and  shape  of  the  active 
plastic  zones,  will  result  as  by-products  of  our  calculations. 

We  develop  a  simple,  approximate  technique  based  on  a  direct  use  of  a  variational 
statement  of  compatibility.  Thus  we  propose  a  trial  function  for  the  stress  function  of  the 
problem,  that  makes  use  of  the  known  asymptotic  information  in  the  near-tip  and  far-field 
limits,  and  depends  cm  arbitrary  parameters  that  measure  the  intensity  of  the  near-tip  fields 
and  other  global  properties  of  the  solution.  Application  of  the  variational  statement  then 
yields  optimal  values  for  these  parameters,  and  hence  determines  die  plastic  stress  intensity 
factor  of  die  near-tip  fields.  The  results  obtained  by  this  novel  method  will  be  compared  to 
the  finite  element  method  results  of  Dean  and  Hutchinson  (1980),  Dean  (1983)  and  Mataga 
(1986)  for  the  linear-hardening  material. 
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2.  Anti-plane  strain  (mode  III) 


2.1  The  plastic  stress  intensity  factor  in  small-scale  yielding 


Figure  1 — Steadily-propagating  aack  in  amall-acale  yielding. 

Consider  the  steadily  propagating  semi-infinite  crack  depicted  in  Figure  1.  Under 
the  small-scale  yielding  assumption,  the  elastic  singularity  field  is  imposed  in  the  remote 
limit  so  that 


x  (r,0)  -  Kel  r*(0)  (2tct)-i/2 


as  r -» oo,  |6|  <  it  (2.1) 


v  (r,6)  — 2  (V  /  G)  Kd  v-(0)  (2nr)~1/2  ■  v- 
where  T~,  V  refer  to  the  well-known  elastic  singularity  solution,  Kel  is  the  elastic  stress 
intensity  factor,  G  is  the  modulus  of  elasticity  in  shear,  V  is  the  propagating  speed  of  the 
crack,  and  x~(8),  v“(0)  are  the  angular  variations  of  the  elastic  fields  in  mode  m, 
normalized  such  that  tj“(0)  ■  1  for  consistency  with  the  definition  of  KeI . 
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On  die  other  hand,  we  require  that  in  the  near-tip  limit 
t  (r,0)  -  KPl  ?»(8)  (2jct)«  ■  t° 

as  r— >0  (2.2) 

v  (r,0)  -  KPl  (V/G)  v°(0)  (2nr)*  /  s  *  v° 

where  x°,  v°  refer  to  the  near-tip  solutions  of  Ponte  Castaiieda  (1985c),  KP1  is  defined  to 
be  die  plastic  stress  intensity  factor  in  small-scale  yielding,  s  *  s(a)  is  the  strength  of  the 
plastic  singularity,  and  we  have  normalized  the  angular  variations  of  these  near-dp  fields, 
which  are  also  functions  of  a,  such  that  t2°(0) » 1. 

It  is  convenient  to  work  with  dimensionless  variables,  and  in  anti-plane  strain  we 
select  the  following: 

I«(VKd)2r 

(2.3) 

S  «(G /tb)(v/V) 
where  is  the  yield  stress  in  shear. 

In  terms  of  these  new  variables  the  remote  fields  can  be  rewritten  as 
r*(r,0)  =  r*(0)  (2rci)~1/2 

(24) 

3T(r,0)  =  -2  v**(0)  (2th)-1/2 
and  the  near-tip  fields  as 

l°(r,0)  =  k  T®(0)  (2;tt)s 

(2-5) 

y?(r,0)  =  k  v°(0)  (2jci)s  /  s 

where  k  is  the  nondimensional  plastic  stress  intensity  factor,  and  is  a  function  of  a  only. 

It  follows  from  this  nondimensionalization  that  the  plastic  stress  intensity  factor 
may  be  expressed  as 

KP>  ■  lim  (2iti)"*  t2(r,0)  -  K(a)  (x0  )1+2s  (Kd)_2*  (2.6) 

r-*0 

and  it  now  remains  to  determine  the  functional  dependence  of  k  on  a. 
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Note  however  that  in  the  elastic  limit  we  must  have  KP1  =  Kel,  which  necessarily 
implies  that 

K(l)«l 

Alternatively,  in  the  perfectly-plasdc  limit  we  expect  KP1  *  and  therefore  that 
K(0)  =  1 

2.2  The  variational  statement  of  the  problem 

CT 

=0 

Figure  2 — Domain  for  the  general  boundary  value  problem. 

We  wish  to  develop  a  simple,  approximate  technique  to  determine  the  functional 
dependence  of  k  on  a.  In  this  connection  we  propose  to  make  use  of  the  following 
variational  statement  of  compatibility: 

Let  A  be  the  domain  of  the  boundary  value  problem,  with 
boundaries  Cj  and  Cy,  shown  in  Figure  2. 

If  the  stress  field  X  satisfies  equilibrium 

V-t-0  in  A  (2.7a) 

and  the  boundary  condition 

ta  *  T  •  n  *  0  on  Cx  (2.7b) 

and  if  the  functional  relation 

JA  G  •  5t  dA  *  Jc»  v*  St,,  ds  (2.7c) 


between  the  engineering  strain-rate  G,  and  the  velocity  field,  prescribed  on 
the  boundary  Cv,  v*.  holds  for  arbitrary  variations  of  the  stress  field  5x, 
that  satisfy  conditions  (2.7a)  and  (2.7b). 

Then  the  strain-rate  field  G  is  compatible 

G  -  Vv  (2.7d) 

and  die  associated  velocity  field  v  meets  the  boundary  condition 

v  a  v*  on  Cy  (2.7c) 


The  constitutive  relation  to  be  used  in  conjunction  with  this  statement  is 

G*I  +GP1  (2.8) 

where  £L  is  the  dimensionless  engineering  strain-rate  vector,  X  is  the  dimensionless 
stress-rate  vector,  and 

0  in  Ae 

fiP1-  (2.9) 

.  (cri-Dde'/ie)!  in  Ap 

is  the  dimensionless  plastic  strain-rate  vector,  which  depends  on  whether  the  material 
particle  is  in  an  active  plastic  zone  Ap,  or  in  a  zone  of  current  elastic  behavior  Ae. 

The  advantage  of  this  particular  choice  for  a  variational  statement  (compatibility 
over  equilibrium)  is  that  with  this  choice  it  is  easier  to  compute  all  the  quantities  involved  in 
the  variational  statement  In  fact,  with  the  introduction  of  a  stress  function  /,  that  ensures 
satisfaction  of  (2.7a)  and  (2.7b),  we  can  easily  compute  all  the  quantities  involved  in  the 
variational  statement  via  simple  differentiation 
I  =  Yx(/e3) 

(2.10) 

I--*l 

Our  approximate  technique  then  consists  in  postulating  the  following  trial  function 


for  / 


/(L8)  -  pd)/°(L0)  +  q(l)  /**(l0) 


(2.11) 


/°(l»0)  -  K«X)  g°(0)  I  (2m)*,  g°(0)  =  -  te°(0)  /  (1+s) 


is  the  stress  function  associated  with  the  near-tip  fields, 

/-(L0)  *  g~(0)  I  (2*I)-1'2.  g“(0)  =  -2  cos  (0/2) 

is  the  stress  function  associated  with  the  remote  fields,  and  where  p(t)  and  q(i)  are 
chosen  so  that  the  resulting  stress  and  deformation  fields  conform  with  the  given 
asymptotic  conditions,  i.e.  equations  (2.1)  and  (2.2).  Hence,  we  must  require 

m  =  1  p(oo)  =  0 

q(0)  =  0  q(°°)  =  1 

Also,  since  we  must  have  that  / = /°  =  /**  for  a  =  1,  we  require  that  p(i)  +  q(i)  =  1  for 
a=  1.  Fen-  simplicity,  we  select  our  trial  functions  p(r)  and  q(r)  such  that  p(x)  +  q(i)  =  1 
for  all  values  of  a.  Thus,  we  take 

p(l)«Y2/(Y2  +  I2)  and  q(l)  =  l2/(Y2  +  I2)  (2.12) 

where  yis  arbitrary. 

Having  specified  /,  we  evaluate  1  and  I,  by  computing  the  appropriate  derivatives 
as  indicated  in  (2.10) 

11  =  /.2=  sin0  ft + 1  -1  cos0  ffi 

12  =  -/j.=-cos0  +  sin0  ffi 

It  * ~fj±=  sin2e  Cfji- 1 -1  ft  - 1  "2  /,ee) "  cos20  (f1  /^e  -I-2  /,e) 
l2=/ji=  cos20  fa  +  sin20  (i -2  /i9e  + 1 -1  /j)  -  sin20  (i"1  -l~2  />0) 

Here  the  computation  of  the  partial  derivatives  with  respect  to  i  is  straightforward,  but 
yields  lengthy  expressions  which  we  will  not  include  here,  and  the  computation  of  the 
partials  with  respect  to  0  involve  the  following  derivatives  of  g°(0) 

(g°)’(6)  -  V»(6)  (g°)"(0)  -  (  xro)’(0) 

which  are  expressed  in  terms  of  quantities  readily  available  from  the  numerical  integration 
of  system  (2.5 )  in  Ponte  Castafieda  (1985c). 


We  then  evaluate  the  variation  of  1  by  assuming  that  k  and  y  are  independent 
parameters.  We  simply  have  that 

5i  =  (fa  /9k)  5«c  +  (fa  /dy)  5y 


ct  c* 


Figure  3— Domain  for  the  variational  lUtemcnt  of  the  iteadily-propagating  crack  in  SSY. 


Next  we  apply  our  variational  statement  to  an  appropriate  domain  A  as  shown  in 
Figure  3.  A  has  been  chosen  for  convenience  to  be  a  semi-circle  of  large  radius  R,  and  is 
composed  of  two  parts:  the  active  plastic  zones  Ap  and  the  zone  of  current  elastic  behavior 
Ae.  Note  that  on  the  crack  face  Cj  T„  =  0  as  required  by  (2.7a),  and  that  on  the  line  ahead 
of  the  crack  Cy 

v*  =  0  =J 


=>  Jcv'S^ds  »0 
On  the  rest  of  the  boundary  CR,  we  assume  that  R  is  large  enough  that  we  can  take 
v*  =  v*.  Then 


v**0(r-1/2)  and  Si,,  =  0(r  *~2) 
Thus  the  variational  statement  reduces  to 


Jcr  v*  8tn  ds  — *  0  as  R-*oo 


JAI-0StfK)dA  +  (a-l-lj/^CleV^I-Ol/aKjdA  =  0 


(2.13) 


JAI*01/9Y)dA  +  (a_1-l) JAp (1*7 Je)  i  •  (5l Idy) dA  -  0 
where  A  is  now  the  upper-half  plane  x2  2  0,  and  Ap  refers  only  to  the  active  plastic  zones. 
Note,  therefore,  that  the  plastic  wake  does  not  enter  directly  in  the  calculation  of  K(a),  and 


also  that  no  integrations  are  needed  beyond  the  quadrature  of  the  variational  integrals.  Note 
further  that  although  the  fields  (and  in  fact  also  some  of  the  integrands  in  (2.13))  are 
singular  and  the  domain  of  integration  is  unbounded,  the  integrals  themselves  are  bounded. 

Equations  (2.13)  amount  to  a  system  of  two  algebraic  equations  in  k  and  y  to  be 
solved  numerically  using  a  secant  method  available  as  a  package  from  I.M.S.L.  The  area 
integrals  were  evaluated  as  iterated  integrals  using  successive  Gaussian  quadratures.  The 
difficulty  in  the  known  singularity  in  r  as  r  -» 0  was  treated  by  subtracting  the  singularity 
from  the  integrand  and  evaluating  its  contribution  analytically.  Alternatively,  the  difficulty 
in  the  known  singularity  in  8  as  0  ->  it,  was  treated  by  a  stretching  technique. 


Figure  4— The  elastic-plastic  boundaries  for  the  steadily-propagating  crack  in  SSY. 

Finally,  let  us  state  for  completeness  the  conditions  that  determine  the  elastic-plastic 
boundaries  shown  in  Figure  4.  The  plastic  loading  boundary  Cq  is  determined  by 

lc2  =  li2  +  I22  =  1 

which  is  interpreted  as  an  algebraic  equation  for  the  radial  coordinate  of  the  boundary  &  as 
a  function  of  6. 

The  elastic  unloading  boundary  Cj  is  determined  by 
lele'  »llll  +  l2l2*0 


which  in  turn  is  interpreted  as  an  algebraic  equation  for  the  angular  coordinate  of  the 
boundary  6  as  a  function  of  i. 


In  anti-plane  strain,  we  will  not  consider  values  of  the  hardening  parameter  for 
which  reloading  takes  place,  and  therefore  we  will  not  discuss  reloading  here,  but  it  is  clear 
that  a  similar  treatment  of  the  plastic  reloading  boundary  C2  and  the  secondary  elastic 
unloading  boundary  C3  (if  necessary)  is  appropriate. 

2.3  Results 

The  results  of  the  numerical  calculations  are  given  in  Table  1,  and  are  depicted  in 
Figures  5  and  6.  It  is  seen  that  for  mode  m,  k(oc)  remains  roughly  constant,  and  y  varies 
dramatically,  especially  near  a  ■  1.  As  we  have  already  remarked,  k(1)  is  expected  to  be 
unity,  but  on  the  other  hand,  y(l)  is  undetermined  by  our  choice  of  p(i)  and  q(i)  (since 
p  +  q  =  1,  independently  of  the  value  of  y),  and  its  steep  behavior  near  a  *  1  is  not  too 
disturbing.  In  the  limit  of  small  a,  the  behavior  of  k  seems  to  be  consistent  with  k  -»  1,  as 
expected,  but  our  results  are  not  conclusive  since  we  do  not  trust  our  trial  function  for  /  for 
very  small  a.  This  is  due  to  the  fact  that  the  limit  as  a  -4  0  of  the  near-tip  Helds  is  not 
uniform,  and  the  range  of  dominance  of  these  Helds  vanishes  in  the  limit  as  a  -*  0  (see 
Ponte  Castaneda,  1985c). 

The  active  plastic  zones  corresponding  to  several  values  of  a  are  shown  in  Figure  7 
using  nondimensional  coordinates  in  the  axes.  It  is  seen  that  the  plastic  zone  gets  shifted 
forward  as  a  decreases.  Also  remark  that  there  is  a  sharp  kink  at  the  intersection  of  the 
plastic  loading  and  elastic  unloading  boundaries,  and  that  the  unloading  boundary  naturally 
intersects  the  plastic  loading  boundary  at  the  place  where  its  slope  vanishes. 

Next  we  will  compare  our  results  to  the  finite  element  method  results  of  Dean  and 
Hutchinson  (1980)  for  linear-hardening  in  mode  III,  as  improved  by  Mataga  (1986). 
Figures  8  are  plots  of  the  logarithm  of  j2  in  the  line  ahead  of  the  crack  versus  the  logarithm 

of  £.  Assuming  that  the  power  singularity  solution  holds  close  to  the  crack  tip  we  have 


I 
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sothat 

lni2slnK*  +  s*lni 

We  fit  a  straight  line  through  those  points  in  Figures  8  that  are  close  to  the  tip  (far  from  the 
elastic-plastic  boundaries  which  show  as  kinks  in  the  graphs),  and  we  compute  the  values 
of  s*  and  k\  Table  2  shows  the  results,  and  how  they  compare  with  our  results  for  s  and 
k.  It  is  found  that  s*  is  not  in  very  good  quantitative  agreement  with  the  theoretically 
correct  s,  particularly  for  the  smaller  values  of  a.  Given  this,  we  can  only  conclude  that  k* 
is  in  adequate  qualitative  agreement  with  tc. 

The  active  plastic  zones  corresponding  to  the  finite  element  results  (crosses)  are 
shown  in  Figures  9,  together  with  our  plastic  zones  (solid  lines)  for  comparison.  It  is 
found  that  our  plastic  zones  have  the  same  general  shape  as  the  finite  dement  plastic  zones. 
However,  it  is  not  clear  which  plastic  zones  are  more  accurate.  The  finite  element  plastic 
zones  probably  suffer  from  inaccuracy  due  to  the  relative  coarseness  of  the  grid, 
particularly  at  the  unloading  boundary  far  away  from  the  tip,  where  the  results  are  very 
sensitive  due  to  the  flat  variation  of  Te  near  unloading.  Alternatively,  our  results  are 
probably  not  very  good  for  the  smaller  values  of  a  due  to  the  expected  decrease  in  the 
range  of  dominance  of  the  near-tip  fields. 

3.  Plane  strain  (mode  I) 

3. 1  The  plastic  stress  intensity  factor  in  small-scale  yielding 

In  plane  strain  the  appropriate  dimensionless  variables  are 
t  ■  (c0  /  KeI  )2  r 

flop  *  Oop  f  °o  C33  *  °33 1  °o 
■  (E  /  C0  )  (va  /  V) 

where  oQ  is  the  yield  stress  in  tension. 


(3.1) 


I 


In  small-scale  yielding  the  remote  fields  are 
fl**(r,0)  =  <r*(6)  (2rci)-1/2 


|0|  <  re 


r*(r,6)  =  -2  v-(0)  (2jc£)-w 

where  <T*(0)  and  V(0)  are  well-known,  and  normalised  such  that  C22**(0)  *  1. 

The  near-tip  fields  are 

fl°(r,0)  =  k  5°(0)  (2jci)* 

(3*3) 

r,0)  *  K  v°(0)  (2jh)!  /  s 

where  s  =  s(a,v),  and  the  angular  variations  have  been  normalized  such  that  O22°(0)  =  1. 
It  is  then  clear  that  for  small-scale  yielding 
k  =  K(a,v) 

and  the  dimensional  plastic  stress  intensity  factor  is  given  by 

KPl  =lim  (2mr*  a^O) »  k(o,v)  (<t0  )1+2»  (K'1)-2*  (3.4) 


Note  that  in  the  elastic  limit  we  must  have 
k(1,v)-1 

and  that  due  to  the  high  level  of  triaxiality  in  the  line  ahead  of  the  propagating  crack  in  the 
perfcctly-plastic  material  (see  Rice,  1980),  we  expect  to  have 
K(0,v)  =  3 

depending  on  the  exact  value  of  v. 


3.2  The  variational  statement  of  the  problem 

In  plane  strain  the  variational  statement  of  compatibility  is: 


and  if 


<*o0,0  -0  in  A 
To«Oo0n0«O  on  Cy 

j*  Da08<7O0dA  »  Jc»v0*8Tads 


(3.5a) 

(3.5b) 


(3.5c) 
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for  arbitrary  So^j  satisfying  (3.5a)  and  (3.5b). 


Then 


Dap  =  (l/2)(va>?  +  vPi0)  (3.5d) 


va ■  va*  on  Cy  (3.5e) 

where  A,  Cj  and  Cv  have  the  same  meaning  as  before. 

The  constitutive  relation  in  dimensionless  form  states 

Dap  *  (1+V)  £ap  -  V  (Sao  +  £33)  +  Dappl  (3 

where  H  is  the  strain-rate  tensor,  £  is  the  stress-rate  tensor,  and 

0  in 

IW1  *  (3 

.  (3/2)  (0-1-1)  (OeV  lie)  Sag  in  A, 

refers  to  the  plastic  strain-rate  tensor.  Here 

fie2  ■  fill2  +  G222  +  fl332  +  3  H\22  “  fill  fi22  “  C22  fl33  “  fl33  fill 
We  can  automatically  satisfy  (3.5a)  by  means  of  a  stress  function  /  such  that 

2n  ■  f 32.  222  *  fT\\  tti2““/o2  (3 

In  addition,  we  require  Q33,  which  is  determined  by  the  plane-strain  condition 


D33  *°! 


£33  -  v  2aa  in  Ac 


£33  ~  v  £oa  +  (3/2)  (or1-!)  (fleV  s«)  £33  in  Ap 


but  this  in  general  involves  solving  a  complicated  differential  equation  in  Ap.  The  difficulty 
is  avoided  by  restricting  our  consideration  to  the  case  of  v  *  1/2,  when 

£33  *  (1/2)  Saa  “  A  (310) 

For  this  reason  we  will  consider  here  only  the  case  of  v  -  1/2.  The  procedure  for 


determining  k  is  then  analogous  to  that  of  the  mode  in  case.  We  take 
/(l8)  -  P(I )/°(L0)  +  q(l)  /**(l>0) 


(3.11) 


5«5CS 
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where 

/°(L0)  -  K  g°(0)  £*  (2*l)»,  g°(0)  .  0*0(0)  /  [(l+s)(2+s)] 

/“(L0)  -  g~(6)  I2  (2jci)-1/2,  g-(0) .  (1/3)  [3  cos  (0/2)  +  cos  (30/2)] 

and  where  p(i)  and  q(t)  are  again  given  by  (2.12). 

Given  /  we  use  (3.8)  to  compute  the  stress  components.  We  also  compute  the 
stress-rate  components  from 

£ii =  ~f.22\  £22  =  ~/.m  £12  =  f.m  £33“  (1/2)  Has 

In  these  calculations  the  following  derivatives  of  g°(0)  are  needed 
(g°)'(0)-- ort°(0)/(l+s) 

(g°)"(0)-  0^(0)  -  (2+s)  g°(0) 

(g°)'"(0)  -  (  o„oy(9)  -  (2+s)  (g°)*(0) 

all  of  which  are  expressed  in  terms  of  quantities  readily  available  from  the  numerical 
integration  of  system  (3.5 )  in  Ponte  Cas  tailed*  (1985c). 

The  variation  of  the  stresses  is  obtained  by  assuming  k  and  y  to  be  independent 
5flop  *  (BOoj/Sk)  5k  +  (dfloj/dy)  &y 
and  again  it  is  found  that  if  we  take  A  to  be  the  upper-half  plane 

Jc,v«*8T«ds  -0 

Thus  the  variational  statement  reduces  to 


Ja  [(^VlJ^aOag-va^  +  WSflaJdA  +... 


(3.12) 

...+  (a-l-1)  JA,  0,7 1.)  [(3/2)  tty  Sflop-  (1/2)  (2aa+  ffssl&IaJ  dA  -  0 
where  Ap  refers  to  die  active  plastic  zones  (both  plastic  loading  and  reloading). 

Again,  for  the  case  of  v  «  1/2,  no  integrations  are  needed  beyond  those  explicitly 


shown  in  the  variational  statement,  which  provides  a  system  of  two  nonlinear  algebraic 
equations  k  and  y  to  be  solved  for  given  a. 
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3.3  Results 

The  results  of  the  numerical  calculations  for  v  «  1/2  are  presented  in  Table  3,  and 
are  shown  graphically  in  Figures  5  and  6.  It  is  observed  that  k  increases  from  unity  as  a 
decreases,  whereas  y  remains  roughly  constant  for  all  a.  More  specifically,  it  is  found  that 
even  though  the  resulting  values  of  k  for  a  «  0.05  and  0.1  seem  to  be  inconsistent  with  the 
rest  of  the  values  of  k,  the  general  trend  of  k  as  a  approaches  zero,  outlined  in  die  figure, 
is  consistent  with  the  perfectly-plastic  result  of  k  *  2.9476  (Slepyan,  1974).  The 
anomalous  behavior  of  die  results  for  k  for  small  a  is  seen  to  be  a  consequence  of  the 
deficiency  in  the  values  of  s  for  small  a  (see  Ponte  Castaneda,  1985c). 

The  active  plastic  zones  corresponding  to  several  values  of  a  are  shown  in  Figure 
10,  where  the  axis  coordinates  have  been  appropriately  nondimensionalised.  In  spite  of  our 
comments  in  the  previous  paragraph,  we  have  included  in  this  figure  the  plastic  zones 
corresponding  to  values  of  a  equal  to  0.05  and  0.1,  for  the  purpose  of  explicidy  showing 
some  plastic  reloading  zones,  which  are  an  important  feature  of  the  plane-strain  mode  I 
problem.  Another  interesting  feature  is  the  fact  that  the  extent  of  the  plastic  zones  in  the  line 
ahead  of  the  crack,  rp,  is  small,  and  in  fact  vanishes  for  a  «1  and  again  fora  -  0.1.  The 
fact  that  rp  vanishes  for  a  *  1  is  intrinsic  to  die  V  *  1/2  case,  but  the  fact  that  it  vanishes 
again  for  a -0.1,  depending  on\the  specific  value  of  v,  is  more  interesting,  and  its  source 
can  be  traced  to  a  change  in  sign  of  the  quantity  jCkr’foO)  -  cOe°(r,0),  which  forces 
j^°(r,0)  to  vanish  when  a  takes  on  the  critical  value. 

Unfortunately,  no  finite  dement  results  are  available  for  die  linear-hardening  plane- 
strain  mode  I  problem  to  compare  with  our  results.  Dean  and  Hutchinson  (1980)  only 
considered  power  hardening  in  mode  I,  and  their  results  are  quite  different  from  ours. 


4.  Plane-stress  (mode  I) 


4.1  The  plastic  stress  intensity  factor 

In  plane  stress  the  appropriate  nondimensional  variables  are  the  same  as  in  plane 
strain  as  given  by  (3.1)  with  the  exception  of  S33,  which  is  set  equal  to  zero.  The  far  fields 
are  as  given  by  (3.2)  where  the  0og~(6)  are  the  same  as  in  plane  strain,  but  the  va~(0)  are 

different  The  near-tip  Helds  have  the  same  form  as  (3.3),  but  the  strength  of  the 
singularity,  and  the  angular  variation  of  the  fields  are  as  given  in  Section  4  of  Ponte 
CastaAeda  (1985c).  The  plastic  stress  intensity  factor  has  then  the  same  general  form  as 
(3.4),  where  again  we  expect  K  =  k(oc,v). 

4.2  The  variational  statement  of  the  problem 

The  same  development  as  for  plane  strain  applies  for  plane  stress,  provided  that  we 
set  Q33  equal  to  zero,  and  we  evaluate  the  pertinent  quantities  that  go  into  the  variational 
statement  (3.12)  using  the  appropriate  near-tip  and  remote  Helds.  Note  that  because  Q33  is 
equal  to  zero,  independently  of  the  value  of  v,  the  variational  statement  applies  equally  well 
for  any  value  of  v.  For  definiteness,  however,  we  will  perform  the  computations  here  for 
v  =  1/2  only. 

Note  that  unlike  s,  which  was  found  to  be  independent  of  v,  k  appears  to  be  a 
functions  of  v  as  well  as  a.  This  is  analogous  to  what  happens  in  the  stationary  plastic 
crack  problem.  Nevertheless,  we  expect  k  to  become  independent  of  v  in  the  limit  as  a 
approaches  zero,  since  it  can  be  shown  (Rice,  1982)  that  the  stress  Held  of  the  perfectly- 
plastic  problem  in  plane  stress  is  independent  of  v. 

4.3  Results 

The  results  of  the  calculations  corresponding  to  the  value  of  v  =  1/2  are  presented  in 
Table  4,  and  in  Figures  5  and  6.  It  is  seen  that  for  plane-stress,  k  increases  slightly  from 
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unity  as  a  decreases,  and  y  remains  roughly  constant  as  a  changes.  Note  that  in  the  elastic 
limit  k  approaches  unity,  but  more  interestingly,  in  die  perfectly-plastic  limit  k  seems  to  be 
consistent  with  the  result  that  one  would  expect  from  the  solution  of  the  perfectly-plastic 
problem  with  a  centered-fan  sector  centered  about  the  crack  line  (see  Rice,  1982  and  Dean, 
1983)  k*  2/V3  s  1.1547. 

The  active  plastic  zones  obtained  are  shown  in  Figure  1 1.  They  have  a  shape  similar 
to  the  corresponding  anti-plane  strain  zones,  but  they  are  taller  and  skinnier,  at  least  for  the 
given  range  of  a. 

Dean  (1983)  produced  the  corresponding  finite  element  method  results  for  this 
problem.  His  results  show  again  that  the  finite  element  method  is  not  very  good  at 
predicting  the  strength  of  the  plastic  singularity  for  small  values  of  a.  However,  we  find 
that  our  plastic  zones  are  similar  in  shape  and  size  to  his. 


5.  Discussion 

It  is  apparent  that  this  new  approximate  method  to  determine  the  plastic  stress 
intensity  factor  of  a  steadily-moving  crack  in  a  linear-hardening  material  is  at  least  more 
efficient  than  the  corresponding  finite  element  method.  It  also  seems  reasonable  that  it 
should  produce  more  accurate  results  for  the  near-tip  fields  as  long  as  the  range  of  validity 
of  the  imposed  near-tip  singularity  field  is  not  too  small.  This  method,  however,  suffers 
from  some  difficultties  that  we  will  address  next 

The  first  of  these  difficulties,  addressed  here  for  convenience  in  the  context  of  mode 
m,  has  to  do  with  the  fact  that  Jj,  as  calculated  from  our  guess  for  /,  becomes  unbounded 
as  0  -» Jt,  whereas  it  is  known  (see  Slepyan,  1973  for  mode  III,  or  Ponte  Castaneda, 
1983b  for  the  analogous  result  in  mode  I)  that  Xi  in  fact  remains  bounded  as  0  -4  it. 
Indeed,  we  find  by  making  use  cf  the  results  for  the  analytical  near-tip  fields  in  the  elastic 


sector  found  in  Ponte  Castafieda  (1985c),  and  the  given  expression  for  /  that  our 
approximate  has  the  following  behavior 

Ii  -  [C  i  p'fc )  (Jt-8)»  +  D  + ...]  i*"1  as  0  -*  n 

where 

iP,(l)--2l2/(E2  +  y2)2 
On  the  other  hand,  we  have  that  the  exact/  is  given  by 
/< r,6)  »/h(r,e)  +  /P«{x2) 

where  /**  yields  a  bounded  contribution  to  Xi  as  6  -»  x,  but  /P*  yields  no  contribution  at 
all,  so  that  die  exact  Xi  remains  bounded  as  0  ->  x. 

Nonetheless,  this  difficulty  should  not  affect  the  values  of  k  and  y  significantly. 
Recall  that  the  stress  lj  itself  becomes  unbounded  as  0  -» x,  and  therefore  the  fact  that 
Xi  *  -lij^  should  stay  bounded  as  0  — >  x  is  a  very  special  circumstance,  that  does  not 

cany  much  weight  in  the  evaluation  of  the  integrals  in  (2.13).  The  reason  for  this  is  that 
even  though  the  approximate  Ij  becomes  unbounded  like  (x-0)*  as  0  -»  x,  the  coefficient 

of  this  term  vanishes  faster  both  as  t  approaches  zero  and  infinity  than  the  next  order  term 
in  8,  which  behaves  like  D  as  0  — >  x,  making  the  artificial  singularity  in  0  unimportant  in 
relative  terms. 

A  way  to  eliminate  the  difficulty,  however,  would  be  to  start  with  a  guess  for  the 
stress-rate  function 

3(r,0)  -  p (£ )  3°(r,8)  +  qfc )  3~(r,8) 

but  this  would  require  integrating  3  with  respect  to  xx  to  obtain  gj.  It  is  conceivable, 
nonetheless,  that  one  may  be  able  to  chose  p(c )  and  q(c )  such  that  the  integration  is  easy, 
but  we  will  not  pursue  this  idea  here  any  further. 

Another  important  difficulty  is  the  fact  that  the  range  of  dominance  of  the  near-tip 
fields,  which  are  a  key  element  of  this  method,  is  expected  to  shrink  to  zero  as  a-*  0.  The 
difficulty  is  compounded  for  the  plane-strain  mode  I  case  where,  as  we  have  seen,  we  trust 
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our  sna’J-a  results  even  less.  Finally,  we  should  add  that  unproved  results  would  also  be 
expected  if  we  included  the  correct  first  order  perturbation  of  the  remote  fields. 

6.  Concluding  remarks 

In  this  work,  we  have  completed  the  knowledge  of  the  near-tip  stress  and 
deformation  fields  of  a  crack  propagating  steadily  and  quasi-statically  in  a  linear-hardening 
material  under  small-scale  yielding  conditions,  by  determining  the  plastic  stress  intensity 
factor  associated  with  these  fields.  The  form  of  these  fields,  and  the  associated  strength  of 
the  singularity  has  been  given  by  Ponte  Castaneda  (1985c). 

The  procedure,  which  so  far  has  only  been  exploited  in  its  crudest  form,  is  fairly 
successful,  provided  that  the  strain-hardening  is  not  too  small,  but  this  is  not  seen  to  be  a 
major  drawback  because  we  expect  that  the  correct  results  for  small  strain-hardening  would 
not  be  too  different  from  the  well-known  perfecdy-plastic  results.  Nevertheless,  the  method 
should  be  subject  to  overall  improvement  by  use  of  a  more  sophisticated  choice  for  the 
stress  function. 

Finally,  it  should  be  pointed  out  that  although  the  present  work  does  not  directly 
address  the  theoretical  prediction  of  resistance  curves,  it  seems  reasonable  that  it  could  be 
extended  to  do  so. 
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Tables 


K* 


1 


a 

K 

r 

0.9 

1.030 

2.895 

0.8 

1.044 

1.880 

0.7 

1.054 

1.524 

0.6 

1.058 

1.318 

0.5 

1.055 

1.163 

0.4 

1.044 

1.006 

0.3 

1.005 

0.685 

0.2 

0.976 

0.513 

0.1 

0.943 

0.509 

-25- 


a 

s 

s* 

K 

k' 

0.9 

-0.483 

-0.475 

1.030 

1.01 

0.5 

-0.394 

-0.375 

1.055 

1.02 

0.3 

-0.325 

-0.297 

1.005 

1.01 

0.2 

-0.276 

-0.247 

0.976 

0.99 

0.1 

-0.207 

-0.169 

0.943 

0.98 

-27 


a 

K 

y 

0.9 

1.031 

1.000 

0.8 

1.063 

1.000 

0.7 

1.093 

1.001 

0.6 

1.122 

1.001 

0.5 

1.148 

1.000 

0.4 

1.169 

1.000 

0.3 

1.179 

1.000 

0.2 

1.173 

0.999 

0.1 

1.137 

0.997 

Table  4 — The  values  of  Kandy  for  plane  stress  mode  I  (v  >1/2). 
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Figure  5 — Plots  of  K  versus  a  for  anti-plane  strain,  plane  strain  (I)  and  plane  stress  (I).  I 
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Figure  6— Plots  of  y  versus  a  for  anti-plane  strain,  plane  strain  (I)  and  plane  stress  (I). 
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Figure  8c— F.E.M.  results  for  the  stress  ahead  of  the  crack  (a  =  0. 1 ) 
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